hep-th/0507189 
AEI-2005-125 
DESY-05-108 
ITEP-TH-49/05 
UUITP-11/05 
ZMP-HH/05-12 



O 

(N 
C 

00 



> 

OS 



Quantum corrections to spinning strings in AdS§ X S 5 
and Bethe ansatz: a comparative study 

Sakura S chafer- NamekP, Marija Zamaklar' 3 and Konstantin Zarembo'il] 



\Q . a II. Institut fur Theoretische Physik der Universitat Hamburg 



Luruper Chaussee 149, 22761 Hamburg, Germany 



00 | sakura.schafer-nameki@desy.de 

t-H ! 

a Zentrum fur Mathematische Physik, Universitat Hamburg 
^ ■ Bundesstrasse 55, 20146 Hamburg, Germany 

13 Max-Planck-Institut fur Gravitationsphysik, AEI 
Am Muhlenberg 1, 14476 Golm, Germany 

Q-i! marzam@aei .mpg.de 

<D . 

7 Department of Theoretical Physics, Uppsala University 
751 08 Uppsala, Sweden 
^ ■ Konstantin.Zarembo@teorfys.uu.se 

«?; 

Abstract 



We analyze quantum corrections to rigid spinning strings in AdS$ x S* 5 . The one-loop worldsheet 
quantum correction to the string energy is compared to the finite-size correction from the 
quantum string Bethe ansatz. Expanding the summands of the string theory energy shift in 
the parameter 1/J 2 and subsequently resumming them yields a divergent result. However, 
upon zeta-function regularisation these results agree with the Bethe ansatz in the first three 
orders. We also perform an analogous computation in the limit of large winding number, 
which results in a disagreement with the string Bethe ansatz prediction. A similar mismatch 
is observed numerically. We comment on the possible origin of this discrepancy. 
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1 Introduction 



Understanding the quantum spectrum of string theory in AdS^ x S 5 is an important open 
problem. Solving this problem will open up venues for testing the ideas of gauge/string duality 
in the genuine stringy regime. It is becoming more and more clear that progress in quantizing 
strings on AdS$ x S 5 is impossible without serious input from the dual M = 4 supersymmetric 
Yang-Mills theory (SYM). One idea that has proved extremely useful on the gauge theory side 
and could potentially be applied to AdS strings, is to compute the spectrum using a Bethe 
ansatz. The Bethe ansatz is the standard approach to quantize integrable systems [I] and it is 
believed that both planar M = 4 SYM and string theory in AdS§ x S 5 are integrable. 

As was observed first at one loop [21 [3] and then at higher orders in perturbation theory 
[2 El E], the planar dilatation operator of M = 4 SYM can be identified with a Hamiltonian of 
an integrable spin chain3. The integrability on the string theory side arises because the classical 
world-sheet sigma-model admits a Lax representation. For the bosonic reduction this almost 
immediately follows [10] from the integrability of the 0(n) model [llj. The Lax pair for the 
full supersymmetric sigma-model in AdS 5 x S 15 [12] was constructed in [T3] . 

Because the classical equations of motion of the AdS string are integrable, their solutions 
can be parameterized by the spectral data of the Lax operator. By reformulating the standard 
solution of the spectral problem [T4] it was shown in [15] that the spectral density for the string 
moving on the R x S 3 subspace of AdS§ x S 5 satisfies an integral equation that strikingly resem- 
bles the large- volume (thermodynamic) limit of the quantum Bethe equations for the spectrum 
of the dilatation operator in the dual gauge theory. These results were extended to other sectors 
[TB"1 WT\ [TBI [19] and eventually to the most general solution including world-sheet fermions [20J. 
Of course the classical approximation in the sigma-model is accurate only at strong 't Hooft 
coupling {i.e. weak worldsheet coupling). In addition, the Noether charges of the string have 
to be large. In order to quantize the string one needs to "undo" the thermodynamic limit and 
turn the integral equations for the sigma-model into discrete, quantum string Bethe equations. 
Such a discretization was first proposed for the su(2) subsector [21], then for other rank-one 
sectors [22J and subsequently for the complete set of Bethe equations with the psu(2, 2|4) sym- 
metry [23]. The quantum string Bethe equations work remarkably well in several tractable 
limits: they have the right classical limit (by construction), reproduce the leading quantum 
corrections for the BMN states and yield the correct energies of massive states in the strict 
strong-coupling limit. 

There are very few explicit calculations for quantum strings in AdS$ x S 5 . One major 
example is string quantization in the plane-wave limit [24J which leads to a solvable string 
theory [22] and can be understood as quantization around the simplest point-like solution of the 
string spinning on S 5 [26]. The curvature corrections [27J to the string states in this background 
(BMN states) were calculated in [21] . Frolov-Tseytlin solutions [23 [30] generalize this setup to 
macroscopic strings and it is possible to quantize fluctuations around these solutions in some 

1 Although the dilatation operator is not integrable beyond leading order in the 1/N expansion [4], the planar 
integrability is still useful in the study of decays of semiclassical strings [7] and in the computation of three- 
point functions [8] . We should also mention that the classical equations of motion of Af = 4 SYM admit a Lax 
representation [3], but we do not know if this property has anything to do with the quantum integrability of 
the planar dilatation operator. 
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cases [3U E2J [33j [33]. For these solutions, the classical string energies can be compared to 
the anomalous dimensions in the gauge theory (see [30, 35J for review), because the 't Hooft 
coupling A combines with the R-charge J into the BMN coupling 1/J 2 = A/ J 2 , which can 
be small even if the 't Hooft coupling is large, provided that the R-charge is large enough. In 
particular, the string action reduces to the effectve action of the spin chain in the limit of large 
J (36]. Generically, one finds that string theory and SYM agree up to two loops and start to 
disagree at three loops. For the quantum corrections the comparison has only been done at the 
one-loop level [371 138] . It would be interesting to understand what happens at higher orders of 
perturbation theory. 

Our goal is to compare quantum corrections to macroscopic strings with the quantum string 
Bethe ansatz at higher loops [HI [22j [23] . The conjectured quantum string Bethe equations were 
rigorously tested at infinite A, but they can potentially receive 1/v^A corrections |2T] . Com- 
parison of the quantum string Bethe ansatz to the direct quantum string calculation provides 
an explicit check of whether such corrections are present at 0(1/ V\) or not. Furthermore, 
the string Bethe equations are known to exactly reproduce the first two orders of the SYM 
perturbation theory independently of J [3T5], and we can just expand the energies computed 
from them in the 't Hooft coupling to find the two loop anomalous dimensions in SYM. In this 
way we can extend the analysis of [371 EH] to two loops. 

Let us briefly review the classical string configurations that we shall study. The one-loop 
quantum corrections were computed for two classes of string solutions - for circular strings 
rotating in S 5 with two independent angular momenta [211 E2] arid for circular strings spinning 
in AdS 3 and rotating around S 5 [33]. The first case is plagued by instabilities [291 [31] an d for 
this reason we shall concentrate on strings moving in AdS% x S 1 C AdS^ x S 5 [30] (throughout 
the paper, we shall adopt the conventions of [33J). The relevant part of the AdS*, x S 5 metric 
in global coordinates is 

d s 2 = - cosh 2 pdt 2 + dp 2 + sinh 2 pd8 2 + d(f) 2 , (1.1) 

where the first three terms are the metric of AdS% and <fi is the angle of a big circle in S 5 . The 
circular string solution has the following form 



p = const , t = kt, 9 = V k 2 + k 2 r + ka, <p = V k 2 + m 2 r + ma, (1-2) 

where 

r 2 = sinh 2 p = — (1.3) 
Vk + k 

£ = , KS +k, (1.4) 
Vk 2 + k 2 



2k£-k 2 = 2VK 2 + k 2 S + J 2 + m 2 , (1.5) 
kS + mJ = 0. (1.6) 

Global charges of the string (the energy E, the spin S, and the angular momentum J) combine 
with the string tension into the following "dimensionless" ratios, which stay finite in the classical 
(A — )• oo, J — > oo, S — > oo) limit 
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Thus 1/y/X or 1/J can be used interchangeably as the loop counting parameters in the sigma- 
model. In addition, at any given order in 1/J one can further expand in the BMN coupling 
1/J 2 = A/ J 2 . In this way one recovers the two-loop perturbative SYM results. 

In section 2 we review the string theory computation and evaluate the energy shift, at 
leading order in 1/J and at the first three orders in 1/J 2 . Although the exact energy shift 
is finite, individual terms of the 1/J expansion diverge. To render the results finite we use a 
particular prescription, the zeta-function regularization. 

In section 3 we compute the energy shift from the quantum string Bethe ansatz, again 
perturbatively in 1/J. Unlike in the string theory calculation, the 1/J expansion is manifestly 
finite. However, the resulting expressions agree with the zeta-regularized string energy shift at 
third order in perturbation theory. 

In section 4 we calculate the energy shift in the non-perturbative regime {i.e., small J) of 
large winding number. The energy shift is finite on both sides in this case. We find a clear 
discrepancy between the Bethe ansatz and the string calculation. In section 5, we present 
numerical results which support the analytical evidence for the discrepancy. 

We discuss our results in section 6. Various technical details are collected in the appendices. 

2 Quantum corrections in string theory 
2.1 Energy shift 

The semiclassical string quantization of [33J yields the following correction to the classical 
energy (jl.4p 

5E string = 5E {0) + 5E osc . (2.1) 
Here the zero-mode contribution is given by 

6E<® 



2 

The oscillator part has the following form 



]- [av + 2k + 2^J k 2 + {1 + r 2 )k 2 - 8Vc 2 + a 2 ^ . (2.2) 



i 00 / 

SE 0SC =- V Wn 2 + v 2 + 2Vn 2 + k 2 - A^{n + 1 ) 2 + a 2 - 4 v / (n-7) 2 + 

K n=l V 

1 4 \ 

7=1 ' 



a 2 



(2.3) 



where the last term is the contribution of the sl(2)-modes, which are the four solutions of the 
quartic equation 

(w 2 -n 2 ) 2 + 4r 1 Vw 2 -4(l + r 1 2 ) (v^+Pw - kn^ = 0. (2.4) 
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The first line corresponds to the transverse and fermionic modes. The various parameters are 
defined as 

v= \[J 2 



72 — m 2 



a 



K + V 2 



k- — 2m 2 — v 1 



m J 



(2.5) 
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K 1 + 



2k 2 k V« 2 + k 2 

2k 2 (l + r 2 ) 



K 2 — V 2 



k 2 — v 2 — 2k 2 r\ 
2(k 2 + k 2 ) 



The sign factors are determined from 



iU) 



ln -n 2 ) Yl^-uj). 



(2.6) 



It is possible to perform a partial summation of the series (12. 3ft . The series is absolutely 
convergent, because the summand decreases as 1/n 2 at n — > oo. Therefore one can sum each 
frequency separately by regularizing the divergences; one adds and subtracts terms of the form 
cin + C2/n before separating various frequencies. This does not change the result, because each 
partial sum is again absolutely convergent. The basic sum is 



n=l 



\/ (n + 7) 2 + a 2 + yjn — j) 2 + a 2 — 2n 



a 



n 



7 



(2.7) 



where {7} denotes the fractional part of 7 and the function F(f3, a) is defined by the following 
integral representation 



F(/3, a) = v 7 " 2 + /3 2 - (3 2 + a 2 / 

Jo 

Using this result we find 



d£ /2Ji(aO 



cosh (3£ — 1 



(2.8) 



5E C 



1 

K 



2F(0, v) + F(0, k) - 4F({ 7 }, a) - 2v - k - 4 7 2 + 4^7 2 + a 2 



1 00 4 / 

n=l 1=1 \ 



« 2 \ 

ooi n — n 

2n / 



(2.9) 



The last sum can be seen to absolutely converge if we use the asymptotic values of the fre- 
quencies ui^ n from [33J. The asymptotic expansion of F((3,a) in 1/a terminates at the second 
order: 



F(P,a) = -a 2 In 



,C-l/2 



a 



(2.10) 



where C = 0.5772. . . is the Euler constant. The dependence on the fractional part of 7 is 
therefore non-perturbative in 1/a and thus in 1/J. In particular it will not be seen in the 
numerical calculations in sec. 5 which will be done for sufficiently large values of J . 
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2.2 Pert ur bat ive expansion 



It is hard to find a useful integral representation for the st(2) modes because of the sign factors 
in ( 12.91) . In computing the perturbative lj J expansion of the string energy shift we shall follow 
a more straightforward approach of evaluating the sum by first expanding all the frequencies 
in 1/ J and then computing the sum order by order in \ j J . As was already observed in [32] 
this procedure is not so harmless, because the sum is not uniformly convergent and modes with 
n ~ J 2 can give a finite contribution. This is reflected in superficial divergences which arise 
starting from second order in lj J 2 . We shall ignore these problems and will use zeta-function 
regularization to sum the divergent series. This approach might not look well motivated but we 
shall find a surprising agreement of this naive summation prescription with the Bethe ansatz 
to third order in 1 / J 2 , which gives us a hint that this prescription may be the correct way to 
compute the energy correction on the string theory side. 

Using the pertubative expressions for the mode frequencies, which are given in appendix B, 
we can write the pertubative expression for the energy shift 5E in powers of 1/ J 2 



~ 8Ei trin 3 

^string \ 

p=l 



J2p 



(2.11) 



It is given by 

SE string = 



I I 00 

-mik — m) H — 2(k — m)m — n 2 + n^/n 2 + 4m(m — k) 



2 y '2 
1 



(2.12) 



n=l 

' ^ 1 1 Is* I O /w\ ^ A 



-m(k — m)(Ak — 11km + 3m 2 

XI 

+ ~{ - 2 ( k - m)m(Ak 2 - 11km + 3m 2 ) + 2(3k 2 - 10km + hm 2 )n 2 + n 4 } 



n=l 



6E string 



n(—A(k — m)m(5k 2 — 15km + 6m 2 ) + 2(k — 3m)(3k — 2m)n 2 + n 4 ) 

8\/n 2 + 4m(m — k) 

— (k - m)m(8k 4 - 52k 3 m + 89k 2 m 2 - 42km 3 + 5m 4 ) 
16 

OO j 

+ Y Yg^ 2 ^ ~ ™)™(8£; 4 - 52£; 3 m + 89A; 2 m 2 - 42/oti 3 + 5m 4 ) 



(2.13) 



n=l 



-(15A; 4 - 128fc 3 m + 279k 2 m 2 - 202km 3 + Um 4 )n 2 

-(15k 2 - 38km + 19m 2 )n 4 - n 6 } 
1 



16(n 2 + 4m(m — k)) 3 / 2 
x {4(k - m) 2 m 2 (A5k 4 - 324k 3 m + 621k 2 m 2 - 370A;m 3 + 60m 4 )n 
-2(k - m)m{53k 4 - A81k 3 m + 1083A; 2 m 2 - 815A;m 3 + 192m 4 )n 3 
+ (15A; 4 - 218A; 3 m + 603fc 2 m 2 - 556A;m 3 + 164m 4 )n 5 
+ (15A; 2 - 44A;m + 25m 2 )n 7 + n 9 } . 



(2.14) 



We shall compare this expression to the energy shift calculated using the Bethe ansatz in the 
next section. 
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3 Bethe ansatz 



3.1 Classical limit 

Classical solutions for the string moving in AdSs x S 1 are uniquely specified by the spectral data 
of the Lax operator. One can introduce the spectral density p(x) defined on a set of intervals 
Ci = (aj,bi). The spectral density satisfies a singular integral equation [16] 



2Jdy^ = 2«kj-2«(^ + ^), xeC, (3.1) 
J x — y \ x — 1 x + ly 

This can be called the classical Bethe equation, as such type of equations arise in the thermo- 
dynamic limit of quantum Bethe equations. 

In addition, the density obeys a set of normalization conditions 

dx — = _27rm, (3.2) 

x 

dx P^l = 2ir(E-S-J), (3.3) 

x z 

Jdxp(x) = 2n(£ + S - J). (3.4) 

Here 2irm is the total world-sheet momentum which must be quantized because of the periodic 
boundary conditions on the world-sheet coordinates. 

We shall consider the simplest solutions of (13. ip with only one cut C = (a, b) which cor- 
responds to the circular string (II. 2p . There is only one mode number k in this case. This 
simplification is crucial and allows us to rewrite the integral equation (]3.ip as an algebraic 
equation for the resolvant 

G(x) = ! dy . (3.5) 

J x-y 

The normalization conditions for the density ( 13 . 2 [) — ( ]3~4|) become boundary conditions for G(x) 

G(0) = 2vrm, (3.6) 
G"(0) = -2n(£ -S - J), (3.7) 
\im zG(z) = 2rr(£ + S-J). (3.8) 

z— >oo 

Multiplying both sides of (13. ip by p(x) /{z — x) and integrating over x we find 



G 2 (z) _ 2n ( k _ 2 ?Z±l!i) G(z) -2n( G{1) + ^—^ G {-\) ) = 0. (3.9) 

z z — 1 / \ z — 1 z + 1 



The boundary conditions (I3.6p - fl3.8p can be used to eliminate G(±l) from this equation. Ex- 
panding (I3.9P at z = and z = oo we get 

kS + mJ = 0, (3.10) 
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in accord with |40j, and 



{J ± m) G(±l) = -irk(£ + S-J)± nm(k + m). (3.11) 

The condition f 1 3 . 1 1) imposes rationality on the spins and requires the integers k and m to have 
opposite signs. We shall assume for defmiteness that m > and k < 0. 
Plugging (13.1 ip back into (I3.9P we get 

(7> _|_ m \ ATT- 2 
k - 2 \ 2 _ x J G(z) + -^—j [k{£ + S-J)z- m{k + m)\ = 0. (3.12) 

The solution of this quadratic equation is 



G{z) =n\k-2 z2 _ x j + z2 _ x , (3.13) 

where 

P(z) = k 2 z 4 - Ak{£ + S)z 3 + 2{2J 2 + 2m 2 - k 2 )z 2 + Ak{£ - S)z + k 2 . (3.14) 
The resolvant determines the density through the discontinuity on the cut 

G(x + i0) -G(x-iO) = 2nip(x), x E C, (3.15) 

and we find 



P(x) = v 2 V • (3.16) 



-P(x) 
x 2 — 1 

We need one extra condition to express the energy in terms of the spin and the angular 
momentum. This condition cannot arise from equation (I3.9p . Instead one should look more 
closely at the structure of the density p{x). For general values of the energy, the angular 
momentum and the spin, the density is real on two cuts, whereas we have assumed that the 
solution has only one cut. This can be made consistent by requiring that the discriminant of 
the quartic polynomial f !3.14j) is zero, then P(z) has one double root (fig.l) 

P(c) = 0, P'( c ) = 0. (3.17) 

These two equations determine the dependence of the energy on the angular momenta, £ = 
£(S, <J), in a parametric form and are equivalent to fll.4j> . ( 11. 5ft upon the identification 

' ' <) . (3.18) 



2 Vc 



3.2 Quantum corrections 

If the integral equation (13. ip is interpreted as the classical limit of some Bethe equations^], the 
density p(x) has the meaning of an asymptotic distribution of Bethe roots in the limit when 

2 Bethe ansatz only works for integrable systems, so here we must assume quantum intergrability of the 
world-sheet sigma-model. There are indeed some indications that integrability is not destroyed by quantum 
corrections |41) . 
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Figure 1: Graph of the quartic polynomial —P(z) (the ordering of the zeroes is a < b < c). 



their number (naturally identified with the spin S of the quantum string state) becomes infinite 

4-7T k — 

^) = ^tE^tt^-^)- ( 3 - 19 ) 

VA k=1 x k i 

The normalization factor 2tc / y/X is the coupling constant of the world-sheet sigma- model. The 
classical (weak-coupling) limit corresponds to A — > oo. Because S scales with y/X according to 
(11.71) . the classical limit coincides with the thermodynamic limit, in which the number of roots 
becomes infinite. 

Our starting point are the quantum Bethe equations proposed in 

/l-^ 1- 



k 7 j^k k 3 x+ Xj V x+x+ x k Xj 



wher^l 



and 



1 

Uk = x k H (3.21) 

x k 

t 1 2wi 
x± + — = u k ±^=. (3.22) 
x k V A 

These equations reduce to (13. ip in the thermodynamic limit when vA, J, S — > oo. Our goal 
will be to compute the leading-order quantum correction to the classical Bethe equations. 

It might seem that ( j3. 201) can only give rise to even powers of 1/ vA, since the equations 
are invariant under y/\ — > —y/X. Nevertheless the odd powers of i/Vx arise in the expansion 
and the leading quantum correction is 0(1/ y/X) for the following reason. The Bethe roots Xk 



3 Although the quantum string can fluctuate in all directions in AdS^ x S 5 , the quantum string Bethe 
equations have the same number of degrees of freedom as in the pure sl(2) sector. On the gauge theory side 
different sectors do not talk to each other because operators with different quantum numbers do not mix [42], 
but it is not a priori clear why various sectors can be separated on the string theory side (see |43j for a more 
detailed discussion of this issue). 

4 Our notation differs from that of [53] by a rescaling of Xk and u^: Xk — > Xk.y/~\/4:Tr, u k — > u k 
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condense into cuts in the thermodynamic limit such that the distance between nearby roots 
goes to zero. But the simultaneous limit of A — > oo and Xk+i — Xk — > is singular in the Bethe 
equations and this singularity gives rise to a local anomaly [44] . The anomaly cancels at the 
leading order j4S], but contributes to the l/y/X quantum correction (37J EE]- We shall calculate 
the anomaly directly from the Bethe equations (13.201) . The calculations are rather complicated 
and the details are given in appendix A. The resulting equation for the resolvant differs from 
(I3.12p by a correction term 



J z + m\ 4tt 



G\z) - 2?r [k - 2 9 - G(z) + 



2 



Z 2 - 1 J w Z 2 - 1 

.2 



k(S + S-J)z-m{k + m) 



4tt z 2 f ^ p'(x)7rp(x) coth7rp(a:) = Q ^ 



y/\ Z 2 — 1 J Z — X 

Solving this quadratic equation we find a density which is of the form (I3.16p . where the 
function P(z) obtains a correction 

5P(z) = f dx Ax)np(x)cothnp(x) (3 24) 

V A n J z — x 

The energy can be found as before, from the requirement that there is only one cut present 

P(c + 5c) + 5P(c + 5c) = 0, P'(c + 5c) + 5P'(c + 5c) = 0. (3.25) 
Expanding the first equation to linear order we get 

?m S£+ ?m Sc+SP{c)=0 . (3.26) 

at ac 

Taking into account that dP(c)/dc = we find 

5P(c) 

58 = -WWd£ ■ (3 ' 27) 



For dP/dE we get from ( 13. 14|) 

9P{C) -Akc(c 2 - 1). (3.28) 



d£ 

Rescaling back to the physical energy we obtain 

irk J x — c 

We can also introduce 

p(x) = - / dff cothf . (3.30) 
K Jo 

Then integration by parts in (13.291) yields 

^Bethe C [ j„ 



§E Betne = / dx J Q 

Trk J [x — c) z 
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Let us see how the one- loop SYM result [37J [38] is recovered. From (I3.17p . (I3.14p we find 
that c = —k/{2J) at large J. Inserting this into ( 13 . 3 1 [) and rescaling x — > AirJ'x, we get for 
the energy shift at the leading order in 1/ J 

in agreement with [37] . 

To perturbatively evaluate the integral (13. 29 p . we shall need to expand various parameters 
characterizing the classical string configuration in a power series in 1/ J . In particular, we need 
to find the zeroes of the quartic polynomial P(x). Recall that P(x) defined in (I3.14p can be 
factorized as 

P{x) = (x - a)(x - b)(x - c) 2 , (3.33) 

For our sign choice (m > 0, k < 0), the roots are ordered as a < b < c. 

The zeroes a, b, c admit an expansion in -j. Solving ( I3.17P perturbatively in 1/J we get 



k k 
2J + 8J r 



-77^ + ^^(2m 2 - Amk + k 2 ) 

1 



+— — (-3m 4 + 16m 3 fc - 23m 2 k 2 + 10mk 3 - k 4 ) + O I ^- } , (3.34) 
8 = (1 — — ) J H — m(m — k) — - m(m — k)(m 2 — 3mk + k 2 ) 



+ 16^ m(m - £;)(m 4 - 7m 3 k + 13m 2 A; 2 - 7mA; 3 + k 4 ) + O [jp j ■ ( 3 - 35 ) 

The expression (I3.35P agrees with the perturbative expansion of the classical string energy 
computed in 



3.3 Mode expansion 

Our starting point is ( I3.29p . which can be written as a contour integral, because the integrand 
has a square-root branch cut along the contour of integration. If we introduce the function 



f(*) = ^£&, (3-36) 
z l — 1 



the energy shift becomes 



iBethe 

J" dz /'(*)/(*) COt(7T/(*)) 



SE Bethe = A _ J K JJK J — — — — — , (3.37) 



k J c , 2ni z — c 



where the integration contour C a b encircles the cut clockwise. We can use the following series 
representation for cot7r/(z) 

oo 

cot(0 = - + 2£^- — . (3.38) 

f ^1 cz _ n z^i 

71=1 
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Inserting this into the contour integral we obtain 



k Jc ab (z-c) k ^ J Cab (z - c){f 2 {z) - n 2 ) ' 

The only singularities of the integrands outside the contour of integration are poles and the 
integrals can be calculated by evaluating the residues. The integrand in the first term has poles 
at z = c and z = ±1. The poles of the second term are at z — ±1 and at z — z n , where the 
z n 's are solutions of 

f(z n ) = ±n, nGN. (3.40) 
Squaring this equation we find that z n 's are the roots of the quartic equation 

P(z) = n 2 (z 2 - l) 2 . (3.41) 

It can be shown that the fluctuation energies around the classical solution are determined by 
the same equation, in accord with the general relationship between fluctuations [IS] and finite- 
size corrections for Bethe ansatz [48]. The residues at z = ±1 are rather complicated, but the 
residues at z = z n are easy to evaluate 

R — = 1 fe) • < 3 - 42 > 

The sign e n of the residue is the same as the sign in the equation f(z n ) = ±n and can be 
determined by analyzing (13.411) with the help of (13.331) 

-1 for z G [— oo, a] U [—1, c] U [1, oo] 

•1 for z G [b, -1] U [c, 1] . 1 ' ' 



3.4 Perturbative expansion and comparison to string theory 

We have evaluated the residues in (13.391) perturbatively in 1/J. The calculations are lengthy 
and are given in appendix C. We also checked that the first two orders are reproduced by a 
direct expansion of the integral (13.311) . Unlike the string sum over modes, its Bethe counterpart 
is manifestly finite at each order of the perturbative expansion. This might indicate that our 
method of computing the series over string modes breaks down at two loops (see also the 
discussion in [22] )• However, if we compare the zeta- regularized sum (12.121) . (I2.13P and (I2.14p 
with the Bethe ansatz, we find complete agreement! We checked this up to the third order 

SE Bethe = §E string^ p = 1, 2, 3. (3.44) 

The agreement at the first two orders implies that the string energy shifts agree with the finite- 
size corrections to the anomalous dimensions at two loops in the SYM theory. At three loops, 
the string Bethe ansatz that was our starting point, differs from the gauge Bethe ansatz [47J 
which computes the anomalous dimensions. 

The agreement between the Bethe ansatz and the direct string calculation is rather spectac- 
ular. The initial expressions look too complicated for this to be a pure accident. Nevertheless, 
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the string and the Bethe calculation have a different status. The Bethe ansatz energy shift 
is automatically finite order by order in \ j J . On the string side we encountered divergences 
despite the complete, unexpanded energy shift being finite. No doubt, there should be a better 
way to approach the weak-coupling (large J) limit on the string side. 



4 Limit of large winding number 

Because of the divergences in the naive 1/ J expansion of the string sum, it would be desirable 
to do an independent test which avoids the convergence issues mentioned earlier. One option 
is to evaluate the energy shifts numerically. This is done in the next section. Here we consider 
a particular regime, the limit of large winding number [\k\ ^> 1), in which the energy shifts 
can be calculated analytically. In this limit J E and m stay finite, but the spin goes to zero: 

5 <C 1. The string remains macroscopic in this limit, since it winds the big circle of S 5 , but in 
AdS§ the string shrinks to zero size (cf. (II. 3ft ). We will have to assume that J /\k\ <C 1, which 
means that there is no overlap with the perturbative regime we have discussed so far. In fact, 
the energy shift turns out to depend on \ jj = \f\/J rather than 1/J 2 in the large-fc limit, 
and it is not possible to compare string quantum corrections to perturbative SYM theory in 
this regime. 

The details of the string calculation are given in appendix D. The result is 



2F (0, V J 2 - m 2 ) + 2F (0, J + m) - AF f{f }, y/J(J + m] 

5E = ^ 

J + m 

VJ + m 



+VmJ + ( J + m) In -j= — -m, (4.1) 

VJ + \fm 

where the function F(f3, a) is defined in (12. 8p . A peculiar property of this result is the depen- 
dence on the fractional part of k/2, which means that the large- k limit of the string energy 
shift depends on whether the winding number k is even or odd. This effect probably arises 
because of the fc-dependent field redefinition of the world-sheet fermions which was used to find 
the spectrum of fluctuations [3~T1 1321 133] . This kind of irregularity does not arise in the Bethe 
ansatz, and also in the zeta-regularized large- J expansion. 



5 In the narrow sense, we are just comparing two mathematical expressions - the string one-loop corrections 
(|2.1[) - (|2.3p and the finite-size correction from the Bethe ansatz (|3.29l) . Each is a well-defined function of the 
parameters k, m and J. If the two expressions agree (or disagree), they must agree (disagree) at all values of 
the parameters, in particular if one of the parameters (k in this case) takes its extreme value. From this point 
of view the limit of large k is just a simplifying assumption that allows us to calculate §E strmg and SE Bethe 
explicitly in some corner of the parameter space. On the other hand, not only the classical energy of the string, 
but also the quantum correction to it stays finite in the large-/c limit. This probably means that the limit of 
large winding (or small spin) is well-defined for this type of string solutions and it would be very interesting 
to study this limit further. The winding number in that, more general setting should be much larger than the 
rescaled quantities £ and J , but should be much smaller than \/A (and thus E and J) in order not to interfere 
with the loop expansion of the sigma-model. 
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4.1 Bethe ansatz calculation 

We begin with the classical limit. To take the large-fc limit it is convenient to rewrite (I3.14p in 
the two equivalent forms 

P(x) = k 2 {x 2 - l) 2 - Ak£x(x 2 - 1) + 4mjx(x ± l) 2 + 4(J q= mfx 2 . (4.2) 

The first two terms blow up in the k — > oo limit unless x is close to 1 or —1. The roots of P, 
a, b and c, thus lie in the vicinity of ±1. Changing the variables to 

x = ±l + p (4.3) 

and taking the limit k — > oo, we get 

P(x) = Av 2 - 8£v + A(J ± m) 2 , at x -± ±1. (4.4) 

Thus two of the roots of P(x) lie near 1 and two lie near —1. The double root should lie at 
x ~ 1, from which we find 

g = j + m (4.5) 

and 

C= 1 - 777- ( 4 -6) 

k 



Solving f)4.4p near x = — 1, we find the endpoints of the cut 



r b \_ 1 (VJ±y^) 
I a J 1^1 



2 



(4.7) 



We see that the cut shrinks to a very small size, whereas the density according to (13. 21) - (13. 4p is 
still normalized to 0(1). Thus the density is highly peaked near —1. Indeed, from (13. 16[) and 
ffl~4"|) we find 

\k\ I 

p( x ) = LAJ2(J + m)v-v 2 - (J -mf. (4.8) 
v v 

The integral (I3.3ip can be easily evaluated in the k — > oo limit. Because the density is large, 
cosh£ in (13.301) can be approximated by 1, and thus 

7T 2 

p = — p , at p — > oo. (4.9) 
2 

We thus get from ( 13.311) 

6E Bethe = 1 f dxp 2, x \ U 1Q \ 

8k J 

Using dx = dv/\k\ and the explicit expression (14. 81) for the density, we find 



J + m yfj + 



5E Bethe = VmJ - In ^ 3- v . (4.11) 

2 VJ^ - \/m 



This clearly disagrees with the string theory calculation (I4.ip . in particular the Bethe ansatz 
result has a regular dependence on k. We shall see this discrepancy also in the numerical 
calculations. Let us also note that even though the explicit computation in this section was 
done in the simplifying large k limit, the deviations between the Bethe ansatz and the string 
theory computation are also observed numerically for finite values of the parameter k (see 
figure |3] in the next section). 
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Figure 2: Energy shifts (SE) x J 2 for J = 5. ..50, m = 3 , k = —2, Bethe vs. semi-classical 
string. 



5 Numerical evaluation of energy shifts 

In this section we numerically compare corrections to the energy of the circular string ob- 
tained by the semiclassical quantization ( 12. 3 p and the one deduced from the proposed quantum 
string Bethe equation ( I3.29p . Both evaluations of the sums are done for various values of the 
parameters. 

We first consider the large- J limit. From figure 2 we see that both functions have the same 
leading order behaviour, in agreement with the earlier analytic results. Next, we try to extract 
the coefficients of the 1/ J 2 expansion of the energy shift numerically. In practice, numerically 
computing higher order effects is hard, since it requires a high numerical precision and stability. 

Yet, by using high precision numerical evaluations let us try to extract the first subleading 
{1/ J 2 ) correction from the exact semiclassical expression ( 12. 3p and compare it with the zeta- 
function regularized result (I2.13p . Subtracting the analytic one- loop piece (12. 12j) from the 
numerical expression for the semiclassical energy shift ( 12. 3 p leads to very unstable numerical 
results, given in table 15.11 



m=3.0 , k=-2. 



J 


50 


100 


150 


200 


250 


(5E strin 3 -SE l )x J 2 


1041 


620 


-82 


-1066 


-2329 


J 


300 


350 


400 


450 


500 


^string _ x j2 


-3871 


-5693 


- 7794 


-10174 


- 12831 



This should be compared to the zeta-function regularized two-loop result ( I2.13P for the same 
values of m and k which gives 

5E 2 = 393.375. (5.2) 

The numerical stability is greatly improved, if instead of subtracting the analytic one-loop 
result (I2.12p . we use the asymptotic numerical value for the energy shift (obtained for J = 10 3 ) 

SE^Ztot = -77.781 . (5.3) 
15 
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Figure 3: Energy shifts (SE) x J 2 for J = 3 , m = 2 , —k = (40. ..1). The upper curve is the 
string calculation. The lower curve is the prediction of the Bethe ansatz. 



The results are given in table ( 15. 4p . We see that it is much less fluctuating compared to the 
result in table ( 15. ip . The deviations from the constant value, may be attributed to higher 
orders in 1/ J 2 and insufficient numerical precision. The average value from the table ( 15. 4p is 
different from the regularized two- loop result (15. 2p . but the numerics is rather unstable and we 
cannot draw any definite conclusions at this point because of insufficient numerical accuracy. 

m=3.0 , k=-2, 

J I 50 I 100 I 150 I 200 I 25T 

(§E atnn 9 - SE asymptot ) x J 2 1170 1167 1147 1120 1087 
J 300 ~ 350 ~400 450 ~500~ 

^string _ fiflasymptot} jj j2 | 1Q4 g | 1Q q 4 | 952 | gQg | §35 

We get much better accuracy if we look at a finite value of J and vary k at fixed m and 
J . We shall take J = 3 and m = 2 and vary k from —40 to —1. The results are given in 
figure [3j The upper curve is the semiclassical string computation, the lower curve is computed 
from the Bethe ansatj§. We see that both the semiclassical and the Bethe energy shifts tend 
asymptotically to constant but different values, which are in good numerical agreement with 
the analytic calculations in the previous section. Here our numerical precision is sufficient to 
discriminate the two results. 

6 Conclusions 

We have compared quantum correction to the energy of macroscopic rigid strings in AdS*, x 5 s 
with the finite-size corrections to the quantum string Bethe ansatz. Taken at face value, the 
two results disagree, but an interpretation of this discrepancy is unclear to us. If we do the 
string calculation in a more naive way by first expanding fluctuation frequencies in 1/J and 

6 By that we mean numerical integration in (I3.29|) . Direct numerical solution of the discrete Bethe equations 
with subsequent extrapolation to the thermodynamic limit requires substantially more involved calculations. 



(5.4) 
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then summing the series over string modes, the straightforward zeta-regularized expansion in 
1/J7 2 agrees with the Bethe ansatz to the first three orders. Perhaps the sum over frequencies 
on the string side should be redefined such that it automatically reproduces zeta-regularized 
\j J expansion. The methods used to evaluate related sums in the context of plane- wave string 
theory |4l5] can be helpful to implement such zeta-function prescription. On the other hand the 
sum is finite and well-defined as it stands and there are no apparent regularization ambiguities. 

Another possible explanation of the discrepancy is that the string Bethe equations receive 
non-trivial l/\/A corrections. We cannot discriminate between these two possibilities at present. 
Studying other classes of string solutions will be certainly helpful to resolve this puzzle. We 
should first of all mention stable circular strings on S 5 which were analyzed both in string 
theory [31] and using the Bethe ansatz |50j . 
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Appendix A Calculation of anomaly 

In this appendix the anomaly term is derived from the quantum string Bethe equations (13.201) . 
The following integral representation turns out to be useful 



in X °^' ■ ■ -'>} = i [" d el £ f -Jjp— V - |4) , (a.i; 



.22L 1 5 

f(x 1 ,...,x s ;x 1 ,...,x s ) J f^{\xl 2 -ldxl x k 2 -ldx k 
where / is an arbitrary function and 

„± 1 



x k + — = u k ± is , (A.2) 



x k 



under the integral (x k on the left-hand-side is defined in (13.221) ). This representation singles 
out a particular branch of the logarithm, so when we write the Bethe equations H 3 . 2 f) in 
the logarithmic form, we should introduce an arbitrary phase which parameterizes different 
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branches of the logarithm 



2tt 

7x 



de < J 



X I 



X i 



™+2 _ -i - 2 

X k 1 X fc 



r + E 



i^fc L 



+ 



x j 



x\ 2 - 1 2 - 1 / X 



k X j 



+ 



x t 



1 a; 



+ 2 



1 / X k Xj 



4 



+ 



(x+ 2 - l)(x-x/ - 1) (XT 2 - 1)(X+XJ - 1) 



(x+ 2 - i)(x+xj - i) (xr 2 - i)(xrx+ - l) 



2tt 



it. 



/y* - 1 - 2 / ^yiH - />"» \ 

^k K^j x j ) 



- 1( + - 



(x+ 2 -l)(x+x, -l)(x+xt-l) 2 - l)(x fe x+ - l)(x fc x, -1) 



/-y» < [ /•>■» 1 ry> 



+ 



x i x k 



(x+ 2 - l)(x+x+ - l)(x fe x+ - 1) (Xj 2 - l)(x+ Xj - l){x k x j - 1) 



2nk. (A.3) 



An important property of this terrible-looking equation is the symmetry with respect to e — > —e, 
which means that the direct strong-coupling expansion starts from order 0(1/ A). The only 
source of i/Va corrections is the first sum over j, in which terms with j ~ k become singular 
in the e — > limit. The contribution of these terms is the anomaly. In the remaining terms we 
can take the limit e — > directly 



Xk x 



Anjxk 47r ^ 

'4^ T + 7t^U-^ 2 -i) 



+ 



2n 



+ 



X 



X : 



1 l 'T* - O' 

1 / X fc 



X, 



+ 



x 



X, 



1 / X k Xj 



J k - ~3 

where we have used the equality 



Ink, 



(A.4) 



U k - Uj 



\Xk Xj^iyXfcXj 1) 



XfcX j 



The next step is to multiply both sides of (1A.4|) by l/(z — x k ) and sum over k. Because of 
the anti-symmetry in k and j, in the double sums l/(z — x k ) can be replaced by 



z- x k 



2 V z Xj z Xj q 



X k Xj 



2(z - x k )(z - Xj) 



Now we can disentangle the "normal" contribution of j — k ~ vA from the local "anomalous" 
contribution of j — k <C vA. In the latter case 



x 



Xj + 



47rx 2 n 



v'Ai./ j - l)p(xj) 



(A.5) 
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according to the definition of the density in (13.191) . Also, 



, x] ( Ann , 

x '^ x ' + —i\7M^ ±ie] (A ' 6) 



"3 ) 



and 



2ie 

1 = =f 



Separating the long-distance contributions from the short- distance ones we find, after some 
calculations 



7 7 4- m \ 4-tt 2 

G\z) - 2tt (k - 2 , ^ G(z) + -5 [k(£ -S-J)z- 2m Jz - m(k + m)\ 

z 2 — 1 / z — 1 



4 2 2x 2 /"v^ i=+oo 

"7x ^3tE (x 2_ 1)( ^_ )2 y * E w +£2 = ' ( A - 7 ) 

v 3 J n=—oo \p £ \Xj) 

where 

The asymptotics of (IA.7I) at z — )• oo shows that the condition (I3.10P does not receive quan- 
tum corrections. Performing the summation in the anomaly term and changing the integration 
variable to £ = \/Xpe/2 we finally get 

,2/ \ „ {, „Jz + m\„, S: 4vr 2 



G\z) - 2tt A; - 2 + - T — ■ + <S - J> - m(fc + m)] 

' cfe = 0, (A.9) 



a/X z 2 — 1 J z — X 

where is defined in (I3.30p . The form of the anomaly used in the main text is obtained 
after integrating by parts and taking into account that 

p' = p'lrp coth 7rp . (A. 10) 

Appendix B Details of string theory computation 
B.l Contribution of s((2) modes 

The main difficulty in evaluating the energy from the string theory is the sum over the roots 
of the quartic polynomial (12. 4p 

6E^ =Y,M9<Cf ) )uj I . n . (B.l) 

The quartic equation is equivalently given by 

u A + a 2 w 2 + aiU) + ao = , (B.2) 
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where 

a 2 = -4k 2 - 2n 2 - Ak 2 r\ - An 2 



fll =8Wfc 2 + /t 2 (l + r|) (B.3) 
a = n 4 -4A;V(l+r 2 ). 



In particular, the absence of the cubic term implies X)/=i ^/.n = 0- The roots can be written 

as 

^l/2,n = 2 ("^n ^ 



(B.4) 



^3/4,n = 2 ( — "^n ^ ^n) , 



where 

Rn = VVl - a 2 



- iii --in 2- ^- (R5) 



and |/i is a real root of the discriminant cubic equation 

V 3 - a 2 y 2 - Aa y + 4a 2 a - a\ — . (B.6) 



That is 



where 



Vl = ia 2 + [M + VM 2 + S 3 ) 1/3 + (m - VM 2 + S 3 ) 1/3 , (B.7) 



5 = ^(-12ao - a 2 ,) 
M = ^(27a 2 - 72a a 2 + 2a\) . 



Furthermore, we need to address the issue of the signs in front of the frequencies. If we 
take all square roots with positive sign, it is clear that for a generic n and J there are two 
possibilities for the relative ordering of the frequencies uj 

I : w 4 < u 3 < uj 2 < u>i (B.9) 

II : 0J4 < u 2 < uj 3 < ui . (B.10) 

In order to discriminate these, consider the large J ^> n limit. The asymptotics are oo\ ~ 
— w 4 ~ 2 J and so (u 2 — n 2 ) > 0. Hence, 

signal) = + 1 > signal) = -1 . (B.ll) 

On the other hand, in the same limit we have u 2 ~ — u 3 ~ n/2J and thus (u 2 — n 2 ) < 0, 
wherefore 

signal) = -1 , signal) = +1 . (B.12) 
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Hence, in the large J limit the eigenvalues are ordered as in the first case in (1B.9[) . Note that 



the ordering of co^ as a function of n keeping J fixed does not change, i.e., the roots do not 
"cross" (see figure [5]). 

Using (IB. lip and (IB. 12ft the expression for 8E sl ^ in the large J limit can be simplified to 
5E m = J2(-co 4 + tu 3 - cu 2 + an) = 2 + u 3 ) = Yl D " + Fn ■ ( B ' 13 ) 

n n n 

In summary, to compute 5E sl ^ one only needs to determine the sum over the combination 
D n + F n . 

B.2 Perturbative expansion of modes 

The combination of sl(2) modes, D n + F n , has the following expansion in \ j J 

5E B{{ - 2) x - / 2k(k -m)+n 2 + n^jAm(m - k) + n 2 \ 1 
— = \<\ 2 )T 2 

4m(k — m){5k 2 — 15km + 6m 2 ) + 2(k — 3m)(3k — 2m)n 3 + n 5 
8^n 2 + 4m(m — k) 

+ -{-2k{k - m){k 2 - 11km + 6m 2 ) - 2(3k 2 - 10km + 5m 2 )n 2 - n 4 
8 

+ (J^{2k(k - m)(k 4 - 23k 3 m + 86k 2 m 2 - 71km 3 + 15m 4 ) 

+ (15A; 4 - 128A; 3 m + 279k 2 m 2 - 202A;m 3 + 44m 4 )n 2 + (15A; 2 - 38km + 19m 2 )n 4 + n 6 } 

+ —r-n \ rrr^ri + 4(Jfe - m) 2 m 2 (45A; 4 - 324A; 3 m + 621A; 2 m 2 - 370km 3 + 60m 4 )n 

16(rr + 4m(m — k)) d ' 2 1 

- 2(k - m)m(53k 4 - 481A; 3 m + 1083fc 2 m 2 - 815£;m 3 + 192m 4 )n 3 

+ (15/c 4 - 218fc 3 m + 603A; 2 m 2 - 556km 3 + 164m 4 )n 5 + (15A; 2 - 44A;m + 25m 2 )n 7 + n 9 }^ 
The other terms, i.e., the transverse and fermionic terms, are as follows 

5E m 



J 4 



(B.14) 



8E 

2k 

= Y.Hk-mf-n 2 )j- 2 

n 

+ —((k- m) 2 (k 2 - 42km - 7m 2 ) + 8(3k 2 - 10km + 5m 2 )n 2 + 4n 4 ) — 

+ — ( - (k - m) 2 (k 4 - 232k 3 m + 962k 2 m 2 - 80km 3 - 11m 4 ) 
128 

- 4(15A; 4 - 260A; 3 m + 594A; 2 m 2 - 340A;m 3 + 23m 4 )n 2 - 16(15fc 2 - 38A;m + 19m 2 )n 4 - 16n 6 )-^ 

(B.15) 
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Note the three-loop term, where the expression at order n 2 has a different structure from the 
one in f ]B.14j) . 



Furthermore, expanding the zero mode part of the energy shift (I2.2p in XjJ we obtain 

111 1 

E^ ' = -m(k — m)— ^ (3k — 7m)(k — m)(k + m) 2 — ; 

2 {J 32 j 

+ —(k - m) (15A; 5 - 135£; 4 m + 182k 3 m 2 - 94k 2 m 3 + 171fcm 4 - 11m 5 ) + O ( ) . 
256 V M ; J7 6 V^ 8 / 

(B.16) 

We shall now combine these terms and obtain the energy shifts up to third order in perturbation 
theory. 

B.3 First and Second order 

The first and second order terms in the \j J 2 expansion of the energy shift (12. II) are 



c x - 2(k — m)m — n 2 + n^Jn 2 + 4m(m — k) 
5E 1 = ^L 9 



2 

n 

r r-iosc ^-(— 4m {k — m) (5k 2 — 15km + 6m 2 ) + 2 (k — 3m) (3A; — 2m) n 2 + n 4 ) 

0-C/ 2 = / , = 

' 8y n 2 + 4m(m — A;) 

+ —((A; -m)(k + mf {-3k + 7m) + 4 (3k 2 - 10km + 5m 2 ) n 2 + 2n 4 ) . 
16 

(B.17) 

The large n behaviour of the summand in 5E° SC is 1/n 2 , which ensures that the energy shift at 
first order is finite. In the second order term the summand has asymptotics 

(5E° sc ) n = — —(k - m) 2 (3£; 2 - 14A;m + 19m 2 ) + O ( X) • (B.18) 
16 \n 2 J 

Thus, there is an anomalous pieces, which needs to be regularized. Applying zeta-function 
regularization the regularized energy reads 

(5E° sc ) reg = — (k - m) 2 (3A; 2 - 14£;m + 19m 2 ) 
9 32 



n(—Am (k — m) (5k — 15km + 6m ) + 2 (k — 3m) (3k — 2m) n + n ) 

8\/n 2 + 4m(m — k) 

+ i(-2 (k-m)m (4k 2 - 11km + 3m 2 ) + 2 (3A; 2 - 10A;m + 5m 2 ) n 2 + n 4 ) J 



(B.19) 
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Combining the zero-mode energy shift with the oscillator contribution, we obtain in summary 
that at order 1/J 2 and 1/J A the shift is 



SE string = 1^ _ m ^ + ^ 2(fc - m)m - n 2 + n^/n 2 + Am(m -~kj 
2 2 

n 

^string _ — m (k _ m)(4k 2 — 11km + 3m 2 ) 



n(—Am (k — m) (5k — 15km + 6m ) + 2 (A; — 3m) (3k — 2m) n + n ) 

+ 4m(m — fe) 

+i(-2 (k-m)m (4k 2 - 11km + 3m 2 ) + 2 (3A; 2 - 10A;m + 5m 2 ) n 2 + n 4 ) J 



(B.20) 



B.4 Third order 



Further expanding the string theory result for the contributions of the oscillators to the energy 
up to third order, i.e., order 1/J7" 6 , yields 

SE° SC = V — \15k 6 - 150Pm + 317A; 4 m 2 - 276k 3 m 3 + 265fc 2 m 4 - 182fcm 5 + 11m 6 
3 ^ 128 1 

n 

+ 4(15k 4 + 4k 3 m - 36£; 2 m 2 - 64£;m 3 + 65m 4 )n 2 - 8(15k 2 - 38km + 19m 2 )n 4 - 8n e } 

+ — r-o \ TT^{^{k - m) 2 m 2 U5k A - 32Ak 3 m + 621A; 2 m 2 - 370A;m 3 + 60m 4 )n 

lQ(n 2 + 4m(m - k)) 3 / 2 1 v ; v ; 

- 2(k - m)m(53k 4 - 481k 3 m + 1083A; 2 m 2 - 815A;m 3 + 192m 4 )n 3 

+ (15A; 4 - 218A; 3 m + 603A; 2 m 2 - 556km 3 + 164m 4 )n 5 + (15fc 2 - 44fcm + 25m 2 )n 7 + n 9 } . 

(B.21) 

The sum is again divergent as the large n behaviour of the summand in (IB.21j) is 



(5E° sc ) n = —(k- m) 2 (5k 2 - 18km + 17m 2 )n 2 

3 1 /IX (B.22) 

+ (k - m) 2 (15k 4 - 2A8k 3 m + 76Qk 2 m 2 - 152km 3 + 91m 4 ) + O — 

128 \n 2 

We again apply zeta-function regularization. In the present case, we need to evaluate the 
Riemann zeta function ((s) = Yl^=i l/n s at s = —2, 0. The values can be calculated by writing 
the zeta-function as 

coo = 73^7 E 2^+i u Y k + ' (R23) 

n=0 k=o v ' 

and evaluating the inner sum first. This results for k > 1 in 

C(-*+l) = ~, (B-24) 
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where Bk are the Bernoulli numbers. Now S3 = and therefore only £(0) gives a non- vanishing 
contribution, namely £(0) = —1/2. The regularized contribution from the oscillators to the zero 
modes is thus 

(5E° sc ) reg = —(k - m) 2 (15A; 4 - 248fc 3 m + 766k 2 m 2 - 752km 3 + 91m 4 ) + ^ - - • , (B.25) 
zoo 

n 

where the dots indicate the non-zero mode contributions, with the terms in (IB. 221) subtracted. 

Combining all terms, we arrive at the third order energy shift as computed from the string 
theory side 

5E stnng = „ _ m ) m ( 8fc 4 _ h2 k 3 m + 89k 2 m 2 - 42km 3 + 5m 4 ) 
16 

+ ^ 2( y k ~ m ) m {^ ~ 52k 3 m + 89k 2 m 2 - 42km 3 + 5m 4 ) 

n 

- (15A; 4 - 128A; 3 m + 279k 2 m 2 - 202km 3 + 44m 4 )n 2 - (15A; 2 - 38km + 19m 2 )n 4 - n 6 } 

+ —r-n \ TTW7^\4(k - m) 2 m 2 (45k 4 - 324k 3 m + 621£; 2 m 2 - 370km 3 + 60m 4 )n 

16(m + 4m(m — k)) 6 ' 1 1 

- 2{k - m)m{53k A - 48lk 3 m + 1083A; 2 m 2 - 815A;m 3 + 192m 4 )n 3 

+ (15fc 4 - 218£; 3 m + 603A; 2 m 2 - 556A;m 3 + 164m 4 )n 5 + (15A; 2 - 44A;m + 25m 2 )n 7 + n 9 } . 

(B.26) 

We shall see subsequently, that this regularized energy shift agrees with the prediction from 
the Bethe ansatz. 



Appendix C Details of Bethe ansatz computation 

C.l Zero- modes 

The zero mode integral is 

SE^ = r<[ dzP^-. (C.l) 
k Jc ab {z - c) 

By deforming the contour to infinity, we pick up the residues at z = c and z — ±1. 

Combining these residues and subsequently expanding them in 1/ J by making use of (13.341) . 
yields 

SE^ = -m(k-m) — 
2 v ' J 2 

1 1 

m(k - m)(4k 2 - 11km + 3m 2 )— (C.2) 

8 yj 

1 1(1 

H (k - m)m(8k A - 52£; 3 m + 89£; 2 m 2 - 42fcm 3 + 5m 4 )— + O — 

16 J b \J S 

Comparison to the string theory result, which were computed in the previous section shows 
that up to third order in the If J 2 perturbation expansion, the zero- mode terms (1C.2j) agree 
with the ones of the zeta-function regularized expressions on the string side. 
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C.2 Non-zero modes 

The non-zero mode contributions come from the sum in (I3.38j) and are 



SE - = E SE » jo E r d mi - w ^ (a3) 



n=l n=l v ; w v ; y 

Again, deforming the contour to infinity, we pick up (possibly non-trivial) residues at z = c, 
z = oo, z = ±1 as well as 2 = z n , where z n were defined in (13.401) . 

The residues at z = c and z = 00 vanish. The residue at z = z n was evaluated in (13.421) . 
In order to expand this in 1/ J, one first needs to solve (I3.4ip perturbatively for z n (note that 
there are two roots z n each for positive n and for negative n). 

The expansion of (I3.42p yields up to third order 

Res 2=2n 
1 



- <j^2fc(m — k) — n + ny n 2 + Am{m — k) j 

' " 2k{m - k){k 2 - 11km + 6m 2 ) + 2(3£; 2 - 10km + 5m> 2 + n 4 

■4(k — m)m(5k 2 — 15km + 6m 2 )n + 2(k — 3m) (3k — 2m)n 3 + n 5 \ 1 



■1 



+ — I -9k 6 + 184k 5 m - 848fc 4 m 2 + 1380fc 3 m 3 - 934fc 2 m 4 + 252km 5 - 20m b 



^/n 2 + 4m(m — k) J J 

\ - 9fc 6 + 184fc 5 m - 848/ 
32 ^ 

- 2(15/fc 4 - 128Pm + 279£; 2 m 2 - 202km 3 + Um 4 )n 2 - 2(15/fc 2 - 38£;m + 19m 2 )n 4 - 2n 

+ 7-5 r tttttt I 4(Jfe - m) 2 m 2 (45A; 4 - 324A; 3 m + 621A; 2 m 2 - 370fcm 3 + 60m 4 )n 

(n 2 + 4m(m - k)) 3 / 2 \ 

- 2(k - m)m(53k 4 - 481k 3 m + 1083k 2 m 2 - 815A;m 3 + 192m 4 )n 3 
+ (15/fc 4 - 218A; 3 m + 603fc 2 m 2 - 556km 3 + 164m 4 )n 5 

1 



+ (15A; 2 - 44£;m + 25m 2 )n 7 + n 9 

+ 



1 

J* 



(C.4) 

Finally, there are the residues at z = ±1, which contribute to the n- independent terms of 
the summands 5E^ 



Res 2= i + Res 



z=-l 



( I,- 2 — m2 )~7p — -^{k — m)(k + m)(k 2 — 8km + 3m 2 )— 4 



- — k(-9k 5 + 152k 4 m - 608k 3 m 2 + 816fc 2 m 3 - 410km 4 + 64m 5 ) — + O ( — 

32 J \\J 



(C.5) 



25 



Putting the residues in (1C.4j) and (1C.5|) together we obtain 
5E {n) 

= — / — 2m[m — k) — n 2 + n\Jn 2 + 4m(m — "fc) \ 

+ i| - 2(jfe - m)m(4k 2 - 11km + 3m 2 ) + 2(3£; 2 - 10£;m + 5m 2 )n 2 + n 4 

— 4(fc — m)m(5k 2 — 15km + 6m 2 )n + 2(/c — 3m)(3A; — 2m)n 3 + n 5 \ 1 



+ —{A(k- m)m(8k 4 - 52k 3 m + 89k 2 m 2 - 42km 3 + 5m 4 



^/n 2 + 4m(m — fc) J J7" 4 

J 4(fc - mWRfc 4 - 52fc 3 r7 

32 



- 2(15/fc 4 - 128fc 3 m + 279A; 2 m 2 - 202km 3 + Um 4 )n 2 - 2(15k 2 - 38fcm + 19m 2 )n 4 - 2n 

+ r TTTTIo ( 4 ~ m) 2 m 2 (4hk 4 - 324k 3 m + 621A; 2 m 2 - 370km 3 + 60m 4 )n 

[n 2 + Amym — k)) 6 ' 2 \ 

- 2(k - m)m{53k A - 481k 3 m + 1083k 2 m 2 - 815A;m 3 + 192m 4 )n 3 
+ (15A; 4 - 218A; 3 m + 603A; 2 m 2 - 556km 3 + 164m 4 )n 5 

+ (15k 2 - 44km + 25m 2 )n 7 + n 9 1 ' ' 



i4 



. j 6 

The complete energy shift is then 

oo 

5E = 5E W + <5£ (n) , (C7) 



(C.6) 



71=1 



where the various terms are written out in (10.20 and (1C.6I) . 

In summary, the Bethe result agrees with the string results (IB. 201) . ( 1B.26|) including order 
1/J 6 . 

Appendix D Details of the large k string computa- 
tion 

We evaluate the energy shift SE strin9 in the large k limit, for fixed m and J. Again, the 
problematic part in the computation are the o;-dependent terms, for which we are forced to use 
approximations for finding the roots in different regions of the parameters. 

Note that first expanding the summands in f |2.2 j) and (I2.3P 1/k before summing them yields 
divergent expressions. However, unlike the divergences that occured in the 1/J expansion at 
second and third order in perturbation theory, these divergences cannot be removed, using stan- 
dard regular isat ion procedures such as zeta-function regularisation as they contain logarithmic 
divergences. The origin of this divergence is the irregular dependence on k of the resummed 
expression (12. 9p . 
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Figure 4: Profiles of the summands for k = 400, k = 100, k = 20 and k — 5, respectively, with 
(J = 3, m = 2). 




Figure 5: Profiles of the u frequencies for k = 400, {JJ — 3, m — 2). The plot on the right hand 
side zooms into the plot on the left hand side. 



In order to ascertain what kind of function we are summing, it is useful to numerically plot 
the summands. This is done in figure 3 for various, mainly large, values of k. Solving ( IB. 21) in 
the limit n ~ \k\ — > oo we find, up to 0(1/ k 2 ) corrections 



Wnl,2 



n±2\k\ ± — 
\k\ 



—n — 



mj + 
(J + m)n 



n±2\k\ (J + m) 



2\k\(n 2 -k 2 ) 



n±\k\ 2 
(J + m)\k\ ± ^(J - m) 2 n 2 + Amjk 2 



(D.l) 
(D.2) 



These expressions approximate the frequencies well enough in the entire range of n, except for 
n — \k\ 
find 



0(1), where 1/k corrections to U2 and U3 blow up. Solving (IB. 21) in that region we 

3|fc|, 



-n, 



\k\ ± y/(n- \k\Y + {J + mf. 



(D.3) 
(D.4) 
(D.5) 



Comparing (ID. 51) to ( ID. II) . (ID. 21) we see that u + asymptotes U2 at n ^> and at n 
while w_ asymptotes at n \k\ and W2 at n <C Thus Co>2 and interchange at n = |/c| 
by passing through the singularity. 

Computing the sign factors from ( 12 .6p we get 



n < 


\k\ : 


signCi = 


1. 


signC^ = 


-1, 


signCf 1 = 1, 


signC^ = — 1, 


n > 


\k\ : 


signCi = 


1, 


signC^ = 


1, 


signCp = -1, 


signC^ = — 1, 


n — 


\k\ ~ 1 : 


signC[ n ^ = 


1. 


signC^ = 


-1, 


signC^ = ±1 . 





(D.6) 
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We are now ready to compute the sum over modes. To do that we divide the range of summation 
into three parts 



(I) l<n<|jfc|-s-l 

(II) \k\ - s < n < \k\ + s 

(III) |A;|+s + l<n, 



(D.7) 



where 1 -C s -C \k\. In the regions (I) and (III) the summation of 0(l/k) terms can be replaced 
by an integration over x = n/\k\ 



oo 4 , 

E E ( si 9 nC i n)ijJ n,I 
n=l 1=1 ^ 



K 2 \ 



|fc|-s-l 



|fc| — s — 1 



n 



dx 



2m J + (J + m 



2 J= 4 £ (\k\-n)-2( J+m y J2 - 

' n=l n=l 

(2 - x 2 )(J + m) + xy/(j - m) 2 x 2 + Am J 



1-x 2 



+ V [2 v // 2 + ( l 7 + m) 2 -3/l + 2(J + m) 2 T dx ( - -) 

l= _ s 1 J Jl-s/\k\ \X - I X J 



k 2 - \k\ + (J + m)Vm~J + mj + F(0, J + m) + {J + m) 2 In 

Combining this with the expansion of the zero modes 

lj-m\ 1J 2 + 10 Jm + m 2 



VJ + m 



(D.8) 



6E® 



1 + 2 



2|fc|( l 7 + m) + 



J + m \ \ J + m 
and (12. 9p we obtain 

k 2 - 4 7 2 + mj + 2F (0, VJ 2 - m 2 ) + 2F (0, J + m) - 4F ({M}, y/jtf + rn)} 



(D.9) 



5£ 



+VmJ + (J" + m)ln 



VJ + m 



(D.10) 



Since 7 = |fc|/2 + 0(l/fc), this expression has a finite — > 00 limit, as was observed numerically. 
In order to determine the asymptotic values of the constant one needs the expression for 7 with 
an 0(1 /k) accuracy 



\k\ m(2J + m) „ ( 1 



(D.ll) 



which implies 



<5£ 



2F (0, VJ 2 ~ m 2 ) + 2F (0, J + m) - AF ({f}, 

VJ + m 



m 



+VmJ +(J + m) In 



J + Vm" 
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— m . 



(D.12) 



For large enough a, the function F(a,(3) can be approximated as in (I2.10p . and thus the 
previous sum can be further simplified to 

SE = -^(J + m) \n(J + m) -(J-m) \n(J-m) + 2J\nJ + 

+V^J - (J + m)\n(^/j + y/m) - m +0 . (D.13) 
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